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ABSTRACT 



Context. The equilibrium structure of sunspots depends critically on its magnetic topology and is dominated by magnetic 

forces. Tension force is one component of the Lorentz force which balances the gradient of magnetic pressure in force-free 

configurations. 

Aims. We employ the tension term of the Lorentz force to clarify the structure of sunspot features like penumbral 

filaments, umbral light bridges and outer penumbral fine structures. 

Methods. We compute vertical component of tension term of Lorentz force over two active regions namely NOAA AR 

10933 and NOAA AR 10930 observed on 05 January 2007 and 12 December 2006 respectively. The former is a simple 

while latter is a complex active region with highly sheared polarity inversion line (PIL). The vector magnetograms used 

are obtained from Hinode(SOT/SP). 

Results. We find an inhomogeneous distribution of tension with both positive and negative signs in various features of 

the sunspots. The existence of positive tension at locations of lower field strength and higher inclination is compatible 

with the uncombed model of the penumbral structure. Positive tension is also seen in umbral light bridges which could 

be indication of uncombed structure of the light bridge. Likewise, the upward directed tension associated with bipolar 

regions in the penumbra could be a direct confirmation of the sea serpent model of penumbral structures. Upward 

directed tension at the PIL of AR 10930 seems to be related to flux emergence. The magnitude of the tension force is 

greater than the force of gravity in some places, implying a nearly force-free configuration for these sunspot features. 

Conclusions. From our study, magnetic tension emerges as a useful diagnostic of the local equilibrium of the sunspot 

fine structures. 
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1. Introduction 

The equilibrium structure of sunspots is obviously domi- 
nated by magnetic forces. Early researches on this problem 
dealt mai nly with the probl em of the global equilibrium of 
sunspots (jMever et al.lll977l ). The sunspot was modelled as 
a magnetic flux rope where the lateral force balance was 
envisaged as a pressure balance between the photospheric 
plasma pressure and th e magnetic -| - plasma pressure in- 
side the flux rope (e.e.. IChitrd 119631 ). Since the magnetic 
field is divergence free, the field lines of the flux rope must 
bend back into the photosphere resulting in a "closed" field 
topology. The resulting curvature of the field lines pro- 
duces a magnetic tension that should basic ally have a down- 
ward vertical component. The model of Me ver et al.l (|1977[ ) 
also put constraints on the field line curvature, since the 
equilibrium becomes unstable when the radius of curva- 
ture is shorter than a certain value. With the availabil- 
ity of high resolution magnetograms of sunspots, it has 
become clear that the earlier models of sunspots might 
no longer be adequate to explain the dynamical equilib- 
rium of various fine structures seen in the umbra as well 
as in the penumbra. One simple diagnostic of the verti- 
cal equilibrium is the vertical component of the magnetic 
tension, which can be determined from the lateral gradi- 
ents of the vector magnetic field. Information about the 
vertical component of the magnetic tension was already 



found to be very useful, as e.g. seen in the correlation of 
low tension forc e with large magnetic shear i n early vector 
magnetograms (jVenkatakrishnan et al. 19931) measur ed by 
the MSEC vector magnetograph (jHaevard et aIlll982D . The 
magnetic tension measurements at the polarity inversion 
lines underlying filaments/prominences are like-wise very 
important since the vanishing of magnetic tension at these 
highly sheared locations makes the prominence structure 
extremely vulnerable to dynamical instabilities, via ther- 
mal instabilities (Venkatakrishnan 1990b). 

Modern observations of sunspots have revealed the 
existence of different fine structures. I n the umbra w e 
have the umbral dots and light bridges ISobotkal (|1989[ ): 
Sobotka et al. (1997b, ^). In the penumbra we have spines 
(stronger, more vertical field) wrapping arou nd the in- 
trasp i nes (weaker, more horizontal field) (JLites et al.l 
I1993I : iBorrero et al.l [2008) . A few models have also been 
proposed to explain these structu res (e.g., uncombed 
model of Solanki fc MontavonI (|1993D . the gappy model of 
ISpruit fc Scharmeij (|200'6[ )T 

In this paper, we investigate the height variation of the 
sunspot fine structure using the calculations of the vertical 
component of the magnetic tension force. The expression 
for computing the vertical component of tension force is 
given in Section 2. In Section 3, we describe the data sets 
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used. Section 4 describes the analysis and results. Finally 
in Section 5 we present our conclusions. 

2. A Brief Description of the Magnetic Tension 
Force 

In any plasma with magnetic field B and plasma pressure p, 
the equation for magneto-hydrostatic equilibrium is given 
by (Parker 1979), 



(V X B) X B/An -Vp + pg^O 



(1) 



where p is the plasma density and g is the acceleration due 
to gravity. The first term in Equation 1 is the Lorentz force, 
second term is the force due to plasma pressure and the last 
term is the force on the plasma due to gravity. We can split 
up the Lorentz force (say F) in two terms as, 



(B-V)B V(B-B) 
47r Stt 



(2) 



The first term in this equation is the tension force (say 
T). The second term represents the force due to magnetic 
pressure. The vertical component of the tension term can 
be expanded in terms of the horizontal derivatives of the 
magnetic field as: 



T, 



4^^ " dx 



B, 



dB, 



B, 



,dB^ dB. 



)] 



(3) 



^ dy "^ dx ' dy 

where, the last component is drawn from the condition, 

V • B = (4) 

The utility of the tension force as a diagnostic of dynami- 
cal equilibrium has not found much atte ntion in the liter- 
ature so far except in a restricted sen se (jVenkatakrishnanl 
I iggOalT bl: I Venkatakrishnan et al.lll993l ). We have computed 
tension force using Equation 3 and expressed it in the units 
of dynes/cm^ . 



3. Data Sets Used 

We have used the vector magnetograms of NOAA AR 10933 
observed on 05 January 2007 and NOAA AR 10930 ob- 
served on 12 December 2006. These data sets are obtained 
from the Solar Optical Telesco pe/Spectro-polar imeter 
(SOT/SP: d Tsuneta et al." "^OOS"; Suematsu et_aL| [2005 
llchimoto et al. 2008: Shimizu ct al. 2008) ) onboard Hinode 
(jKosugiet aL[l2007D . 

The Hinode (SOT/SP) data have been calibrated 
by the standard "SP_PREP" routine developed by 
B. Lites and available in the Solar-Soft package. The 
"SP_PREP" determines the thermal shifts in the spectral 
and slit dimensions and also applies the drift correc- 
tions for calibrating the data from levelO to levell. 
The prepared polarization spectra have then been in- 
verted to obtain vector magnetic fie l d components using 
an Unno-Rachkowsky (jUnnol 119561: iRachkowskvl 119671 ) 
inversion under the a ssumption of Milne-Eddington 
(ME) atmospher e ([Landolfi fc Landi Deel'Innocentil 



Il982t [Skumanich &Jited 119871). We have used the 



"STOKESFIT" inversion code B which is available in the 
Solar-Soft package. The latest version of the inversion code 



is used which returns the intrinsic field strengths along with 
the filling factor. The effect of polarimetric noise on the 
estima tion of vector fields is almost negligible (Tiwa ri et al.l 
l2009al ) and (Gosain, Tiwari and Venkatakrishnan, 2010, 
ApJ, (sent)). The 180° azimuthal ambiguity in our data 
sets have been removed by using acute angle metho d 
(|Harvevlll969l: ISakurai et allll985l: iCuperman et al.lll992[ ). 
The data sets used, have spatial sampling of ~ 0.32 
arcsec/pixel observed in the "fast mode" scans of SOT/SP. 
A simple (NOAA AR 10933) and a complex (NOAA AR 
10930) sunspots are analyzed. The noise in the d ata has 
been m inimized in t he simi la r way as done in [Tiwari 
f |2009[) : iTiwari et all (|2009b[l : [Venkatakrishnan fc TiwaTi 

mm. 



4. Data Analysis and Results 



Tension Contours on Continuum l\/lap(NOAA10933:05jan07:0708hrs) 



200- 




50- 



The code has been developed by T. R. Metcalf 



Fig. 1. Contours of tension forces overlaid on the con- 
tinuum map of NOAA AR 10933 (S04E05). Blue (red) 
colors show negative (positive) contours of ±1.2, ±4, ±12 
millidynes/cm^ . The position of a light bridge(LB) and 
three examples of bipolar sea serpent regions are shown 
by an arrow and circles respectively. An arc is shown for 
which the scatter plots of tension, field strength and incli- 
nation are shown in Figure 2. The heliocentric angle is 9 — 
8° and an arrow points towards the center of the solar disk 
(CSD). 



We calculate the vertical component of the tension force 
for all locations in the sunspots. Figure 1 shows a con- 
tinuum picture of the NOAA AR 10933 overlaid by the 
contours of the magnetic tension expressed in units of 
dynes/cm'^. The red contours depict positive values, while 
the blue contours mark negative values. Several interest- 
ing points can be noted. The most obvious result is the 
existence of locations of positive tension, contrary to the 
expectation from simple sunspot models having a "closed" 
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Fig. 2. Upper panel: The scatter plot between magnetic 
tension and field strength along the arc shown in Figure 
1. Lower panel: The scatter plot between magnetic tension 
and field inclination along the arc shown in Figure 1. 



field topology. Of further interest is the fact that these lo- 
cations of positive tension lie along linear structures in the 
penumbra, which resemble the penumbral filaments. The 
magnitude of the tension is very much comparable to the 
force due to gravity and exceeds the solar gravity in some 
places. For a more detailed comparison, we show in Figure 
Al (online), a map of the magnetic field strength overlaid 
with the contours of vertical tension. In Figure A2 (online) , 
we present a map of the inclination overlaid with contours 
of the magnetic tension. Here again, the locations of posi- 
tive tension largely coincide with regions of lower magnetic 
field strength and larger magnetic inclination. In Figure 2, 
we show scatter plots of field strength versus tension as 
well as inclination versus tension. These scatter plots show 
a clear trend of the positive tension being associated with 
weaker and more inclined fields. 

Apart from the penumbral features, we can also see sev- 
eral umbral features including a light bridge at the western 
portion of the umbra. Even in this case, the regions of pos- 
itive tension coincide with locations of weaker and more in- 
clined fields. Finally, we notice several locations of upward 
directed tension at the edges of the filamentary penumbral 
structures. In Figure 3, we show the contours of tension 
overlaid on the continuum map of NOAA AR 10930. Apart 
from the similar structures of positive and negative tensions 
as seen in the NOAA AR 10933 (Figure 1), the polarity in- 
version line of NOAA AR 10930 also shows high magnitudes 
of tension (10^^ dynes/cm'^). 

5. Discussion and Conclusions 

The variation of magnetic field parameters wi t h op- 
tical depth can be obtained ( CoUado s et al.l Il994t 
IWestendorp Plaza et all l2001bllal: IB eck 200j}_using inver- 
sion schemes like SIR (jRuiz Cobo fc del Toro Iniestalll99a . 



199 4 iDel Toro Iniesta fc Ruiz Cobdll996l: Idel Toro Iniestal 
2003). With the availability of high quality stokes pro- 
filcs from the Solar Optical Telescope ( JTsuncta ct al. 200a 
Suematsu erall 120081: Hchimoto et allEoOSl: IShimizu et al.l 



2008D aboard Hinode (jKosugi et al.ll2007() . these inversions 
are now possible with unprecedented quality. In particu- 
lar, this was used to demonstrate the variation of magnetic 
field o f penumbral filamen ts as a function of the optical 
depth (jBorrero et al.ll2008l) , which showed the wrapping of 
the spine around t h e intr aspine. An inspection of Figure 
2 of iBorrero et al.l (|200^ shows that the magnetic field 
lines can be inferred to have an upward directed curva- 
ture above the intraspine. In the case of our vector mag- 
netograms, the inversions provide the magnetic parame- 
ters at a single optical depth in the atmosphere. It has 
been recognized that the simple inversions based on opti- 
cal depth independent magnetic parameters are consistent 
with a weighted r nean of the optical depth depe ndent in- 
versions like SIR (jWestendorp Plaza et al.l[l998[ ). We can 
thus conjecture that our locations of positive tension are 
manifestations of the upward directed field curvatures that 
would be seen in SIR inversions. This conjecture must be 
verified with detailed SIR inversions over the entire sunspot, 
which will be presented in a future detailed paper. For the 
present, we can only note that the locations of positive ten- 
sion in the penumbra appear to be in the regions of lower 
field strength and larger inclination, which are the proper- 
ties of the intraspine. We can thus conclude that the lateral 
variation of the magnetic field, as inferred from simplistic 
inversions of the stokes profiles, appear to show an upward 
directed tension, which may be in fact located above the 
intraspine. 

Likewise, the upward directed tension seen in 
light bridges could well be due to a similar bend- 
i ng of field lines around a more horizontal field 
(jSankarasubramanian fc Hagenaaill2007l ). Finally, the pos- 
itive tension seen at the edges in bipolar sections of the 
filamentary penumbral structures could well be a sig nature 
of the "sea-serpent" (jSainz Dalda fc Bellot Rubiol l2008[) 
where the penumbral fields dip into the photosphere due 
to hydrodynamic forces caused by convection. 

Apart from the sign, the magnitude of the vertical ten- 
sion also has some very interesting implications. As men- 
tioned earlier, the magnitude of tension attains values upto 
10^^ dynes/cm'^ in some places. This value is comparable to 
the force of gravity which ranges from 5x10"'^ dynes/cm^ 
for quiet Sun density to 10~^ dynes/cm'^ for umbral den- 
sities. This must mean that non-magnetic forces alone will 
not be able to balance this tension force. The only other 
force that can match this force is the gradient of mag- 
netic pressure. As an example, we can use the values of 
spine (~ 2000G) and intraspine (~ lOOOG) of lBorrero et al.l 
(2008) and a scale height A of ^100 km for the vertical vari- 
ation of field strength, to estimate a downward magnetic 
pressure gradient (■^^) of ~ 1.2 x 10~^ dynes/cm^. This 
value is of the same order as that of magnetic tension ob- 
tained in our analysis. It is well known that configurations 
where the magnetic tension balances the magnetic pressure 
gradient are the so-called force-free configurations. Thus, 
the use of high resolution magnetograms leads to the impor- 
tant result that the existence of upward magnetic tension in 
sunspot penumbral fine structure could be a manifestation 
of uncombed magnetic field in force-free equilibrium. 
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Since the observations refer to values at a given optical 
depth, the resultant magnetogram is not at a single geomet- 
rical height. What could be the effect of this corrugation in 
the height levels? One can estimate that the relative am- 
plitude of tension fluctuations (AT/T) will be proportional 
to the relative amplitude in the field variations {AB/B). 
Assuming a vertical gradient of magnetic field of IG/km 
in the penumbra and a corrugation amplitude of 100 km 
for the optical depth, we obtain AB of ~ 100 G. Knowing 
that the penumbral fields are ~ 1000 G, we can estimate 
the relative field variation as 10%, which then is also the 
estimate for the relative variation in the tension force. The 
observed spatial variation in tension is much greater than 
this estimated artefact in tension fluctuations. 



Tension Contours on Continuum Map(N0AA1 0930:1 2Dec06:0304lirs) 




Fig. 3. Contours of tension forces overlaid on the contin- 
uum map of NOAA AR 10930 observed on 12 December 
2006 at an heliocentric angle oi 9 = 31°. Blue (red) 
colors show negative (positive) contours of ±1.2, ±4, ±12 
millidynes/cm'^. Polarity inversion line is shown by dashed 
yellow line. An arrow points towards the center of the solar 
disk (CSD). 



We must remember that the 180 degree ambiguity is not 
an issue whenever the transverse magnetic vector makes an 
angle less than 90 degrees with the potential held. In all 
such cases, the ambiguity is properly resolved by the acute 
angle method. The NOAA AR 10933 is away from the po- 
larity i nversion line and is not highly sh eared (JTiwari et al.l 
l2009bt IVenkatakrishnan fc Tiwajlll2009t ). We note that the 
locations of positive tension discussed thus far all have 
transverse vectors which are not perpendicular to the po- 
tential fields. Hence they have no problem as far as the 
acute angle method is concerned. We do see however that 
this will be an issue for highly sheared portions of the mag- 
netic field. We normally expect low values of magnetic ten- 
sion at such highly sheared regions. But in the case of 
AR 10930, there is strong magnetic tension seen in the 



highly sheared regions at the polarity inversion line (Figure 
3). In such cases, we need to have a careful look at the 
180 degree ambiguity resolution before we can decide on 
the sign of magnetic tension. However, the upward tension 
could indee d be a signature o f the forces that drive flux 
emergence (jLites fc Low! 119971 ). This conjecture is consis- 
tent with the copious emergence of flux observed on 12 Dec 
2006 at the p olarity inversion line (PIL) of the AR 10930 
(jSchriiver et al. 2008). In contrast, the PIL of AR 10930 on 
14 December 2006 in Figure A3 (online) shows a reduced 
value of the tension corresponding to a time of negligible 
flux emergence. Thus, the study of magnetic tension in var- 
ious types of sunspot fine structures promises to yield new 
and exciting information on the equilibrium and dynamics 
of these structures. 
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Appendix A: Additional figures 



Tension Contours on Field Strength Map(NOAA10933:05jan07:0708hrs; 



200 -: / 



Tension Contours on Continuum Map(NOAA1 0930:1 4Dec06:0001hrs) 
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Fig. A. 3. Same as that of Figure 5, except that this sunspot 

„. . , „ ;i ; r TT 1 . ,1 , .1 IS obscrved on 14 December 2006 at an hehocentric anaic 

Fig. A.l. Same as that or l^igure 1, except that the f a _ qo 

continuum image is replaced by the field strength map 

of NOAA AR 10933. Contour levels are ±1.2, ±4, ±12 

millidynes/cm'^ . 



Tension Contours on Inclination Map(NOAA10933:05jan07:0708lirs) 
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Fig. A. 2. Same as that of Figure 1, except that the con- 
tinuum image is replaced by the field inclination map of 
NOAA AR 10933. 



